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ABSTRACT 

We exploit recent constraints on the ages and metallicities of early-type galaxies in 
the Sloan Digital Sky Survey (SDSS) to gain new insight into the physical origin 
of two fundamental relations obeyed by these galaxies: the colour-magnitude and the 
Mg2-oy relations. Our sample consists of 26,003 galaxies selected from the SDSS Data 
Release Two (DR2) on the basis of their concentrated light profiles, for which we have 
previously derived median-likelihood estimates of stellar metallicity, light-weighted 
age and stellar mass. Our analysis provides the most unambiguous demonstration 
to date of the fact that both the colour-magnitude and the Mg2-ay relations are 
primarily sequences in stellar mass and that total stellar metallicity, a-elements-to-iron 
abundance ratio and light-weighted age all increase with mass along the two relations. 
For high-mass ellipticals, the dispersion in age is small and consistent with the error. At 
the low-mass end, there is a tail towards younger ages, which dominates the scatter in 
colour and index strength at fixed mass. A small, but detectable, intrinsic scatter in the 
mass-metallicity relation also contributes to the scatter in the two observational scaling 
relations, even at high masses. Our results suggest that the chemical composition of 
an early-type galaxy is more tightly related to its dynamical mass (including stars and 
dark matter) than to its stellar mass. The ratio between stellar mass and dynamical 
mass appears to decrease from the least massive to the most massive galaxies in our 
sample. 
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1 INTRODUCTION 

The observed properties of early-type galaxies obey sev- 
eral fundamental relations, which have long been thought 
to hide important clues about the physical processes 
that influenced the formation and evolution of these 
galaxies. For example, luminosity, central velocity dis- 
persion o v , mean surface brightness, brightness profile, 
colours and Mg2 absorption- line index all appear to be 
tight l y related to each other in early-type galaxies feauml 
19591: iFaberfc JacksorJll976t IVisyanathan fc Sandagdll977t 
Kormcndvj Il977t iBender et alJ [1993; Dio rgovski fc Davis! 
198711 . The colour-magnitude and the Mg2-oy relations are 
particularly interesting in that they connect the luminous 
and dynamical masses of the galaxies with the physical prop- 
erties of their stellar populations. The tightness and homo- 



* E-mail: gallazzi@MPA-Garching.MPG.DE 



geneity of these two relations must be telling us something 
fundamental about the epoch and the process of for mation 
of early-type galaxies (see for a review lRenzinil l2006 1 . 

The physical interpretation of the above observational 
relations is still a subject of debate. The colour-magnitude 
relation is often interpreted as a sequenc e of increas- 
ing metallicity with increasing luminosit y J^berJ |l9z3; 
[ Worthed Il994t iKodama fc Arimotol Il997l iKodama etld] 
1999). However, it has been proposed that, in addi- 
tion to met allicity, age could at least in part drive 
the relation jG onzalez et alJ |1993|: [T erlevich et alJ Il999l : 
iFerreras et ail Il999t iPoggianti et all 1200 lH . Studies of 
evolution with cosmic time indicate that the slope of 
the colo ur-magnitude relatio n has changed little since 
z ~ 1 iKodama fc Arimotol Il997t IKodama et alJ Il998l : 
IStanford et alJll998t iBlakeslee et al.ll2003l1 . This has been 
used as an argument against age as the primary driver of 
the relation, since the colours of young stellar populations 
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evolve faster than those of old stellar populations. Simi- 
lar interpretations have been proposed for the relation be- 
tween Mg2 index strength and central velocity dispersion 
ay, which is generally thought to ari se from a combina- 
tion of age and metallicitv variations ijColless et alJ Il999t 
iTrager et a.l]l2000al: iKuntschner et aljEioih. 

The difficulty of obtaining unambiguous constraints on 
the relative influence of age and metallicity on the colour- 
magnitude and the Mg2-ov relations is a consequence of 
the difficulty of deriving accurate ages and metallicities for 
large samples of early- type galaxies: age and metallicity 
both tend to redden the colours and strengthen the Mg2 
absorption line in similar ways (e.g., IWorthevlll99il . In 
addition, at fixed metallicity, the a-element-to-iron (a/Fe) 
abundance ratio appears to be larger in the most massive 
early-type galaxies than in the near by stars used to cali- 
brate age and meta llicity estimates fWorthev et alJ Il992t 
^azrjgkis^^ilj |22pJ) . This has been explored recently by 
ThrjmasetHuT<l2004ll . who used 'closed-box' chemical evo- 
lution models with variable heavy-element abundance ratios 
to analyse a heterogeneous sample of 124 nearby early-type 
galaxies, in both low-density and high-density environments. 
According to these models, massive galaxies formed earlier 
and more rapidly than low-mass galaxies, while both the 
colour-magnitude and the Mg2-oy re lations are p rimarily 
driven by metallicity. iBernardi et alJ ll2003allblcld|) carried 
out a more observationally oriented analysis on a sample of 
nearly 9000 early-type galaxies from the Sloan Digital Sky 
Survey (SDSS). They showed that the colour-magnitude re- 
lation reflects a dependence of both colour and luminosity 
on velocity dispersion. By assuming that luminosity traces 
metallicity an d that the scatt e r in c olour at fixed luminos- 
ity traces age, IBernardi et alJ J2005f) explored how age and 
metallicity may be related to velocity dispersion in early- 
type galaxies. 

In this paper, we re-examine the physical origin of the 
colour-magnitude and the Mg2-ov relations using a differ- 
ent approach. Our starting point is a set of statistical es- 
timates of light-weighted age, stellar metallicity and stellar 
mass for a large sample of 26,003 early-type galaxies drawn 
from the SDSS Data Release Two (DR2). We derived these 
constraints in earlier work by using a compre hensive library 
of mo del spectra at medium-high resolution l|Qallazzi et alJ 
120051 . hereafter Paper I), to interpret the strengths of five 
spectral absorption features with negligible dependence on 
the a/Fe ratio. We use here this dataset, together with an 
observational tracer of the a/Fe ratio, to demonstrate unam- 
biguously that both the colour-magnitude and the Mg2-ov 
relations of early-type galaxies are primarily sequences in 
stellar mass and that both the total stellar metallicity and 
the a/Fe ratio increase with mass along the two relations. 
Light-weighted age increases from the least massive to the 
most massive early-type galaxies, with a larger spread at low 
masses that dominates the scatter in the colour-magnitude 
and Mg2-ov relations. The small intrinsic scatter in metal- 
licity at fixed mass also contributes to the scatter in the two 
scaling relations. 

We present our sample in Section [2] below, along with 
a brief description of the method adopted in Paper I to de- 
rive statistical estimates of the ages, metallicities and stellar 
masses of the galaxies. The influence of these parameters on 
the colour-magnitude and the Mg2-ov relations is explored 



in Sections 13 . 1 1 and 13.21 respectively, and their possible de- 
pendence on galaxy environment is addressed in Section 13.31 
In Section E31 we discuss some implications on the relations 
between physical parameters and dynamical mass. Section^] 
summarises our conclusions. Throughout the paper we use 
fi m = 0.3, fi A = 0.7 and H = 70 km s" 1 Mpc" 1 



2 OBSERVATIONAL SAMPLE 

We select our sam ple from the main spe ctroscopic sample 
of the SDSS DR2 jAbazaiian et al.ll2004l) . The SDSS is an 
imaging and spectroscopic survey of the high Galactic lat- 
itude sky, which will obtain u, g, r, i and z photometry of 
almost a quarter of the s ky and spectra of at least 700,000 
objects iYork et al.ll2000h . The spectra are taken using 3"- 
diameter fibres, positioned as close as possible to the centres 
of the target galaxies. Stellar metallicity, light-weighted age 
and stellar mass estimates are available from Paper I for 
a sample of 175,128 galaxies with Petrosian r-band magni- 
tudes in the range 14.5 < r < 17.77 (after correction for 
foreground Galactic extinction) and in the redshift range 
0.005 < z ^ 0.22. The lower redshift limit allows us to 
include low-luminosity galaxies (corresponding to a stellar 
mass of ~ 10 8 Mq), while still avoiding redshifts for which 
deviations from the Hubble flow can be substantial. The up- 
per limit corresponds roughly to the redshift at which a typ- 
ical 10 11 Mq galaxy is detected with median S/N per pixel 
greater than 20. 

We select early-type galaxies on the basis of the light 
concentration index C = Rgo / R50 , defined as the ratio of the 
radii enclosing 90 and 50 percent of the total Petrosian r- 
band luminosity of a galaxy. This parameter has been shown 
to correlate well with morp hological type iStrateva et alJ 
l200lt Ishimasa ku et aljfeoOlT) . Thus, it allow s a rough mor- 
pholog ical classification of SDSS galaxies. IStrateva et alJ 
(2001) propose a cut at C = 2.6 to separate early- from 
late-type galaxies. To limit the contamination by disc galax- 
ies with large bulges, we define here as 'early-type' those 
galaxies with concentration index C 2.8. In this way, 
we select 67,411 early- type galaxies in the redshift range 
0.005 < 2 ^ 0.22. We note that we decided against further 
limiting the contamination of our sample by systems with 
residual star formation by imposing a lower cutoff in 4000-A 
break strength or an upper limit on emission lines equiva- 
lent width. Since these quantities correlate with colour, such 
limits would introduce an unwanted cutoff in the colour- 
magnitude relation. 

Bayesian-likelihood estimates of the r-band light- 
weighted ages, stellar metallicities and stellar masses of the 
67,411 galaxies in our sample are available from Paper I. 
These estimates were deriv ed by comparing the spect rum of 
each galaxy to a library of iBruzual fc Charlotl j2003T) mod- 
els at medium-high spectral resolution, encompassing the 
full range of physically plausible star formation histories. In 

1 IBernardi et alJ 120051) adopted a different criterion to select 
early-type galaxies from the SDSS. They considered to be of early 
type those galaxies for which the r-band surface brightness profile 
is better described by a de Vaucouleurs law than by an exponen- 
tial law (photometric parameter fracDev > 0.8) and that do not 
have emission lines (spectroscopic parameter eclass < 0). 
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practice, we compared the strengths of five spectral absorp- 
tion features in the spectrum of each observed galaxy to the 
strengths of these features in every model spectrum (broad- 
ened to the observed velocity dispersion) in the library. We 
used D4000, H/3 and H<5a±H7a as age-sensitive indices and 
[Mg2Fe] and [MgFe]' as metal-sensitive indices, all of which 
depend negligibly on the a/Fe ratio. 2 This comparison al- 
lowed us to construct the probability density functions of 
age, metallicity and stellar mass for every galaxy. The esti- 
mate of each parameter is given by the median of the cor- 
responding probability distribution, while the Ha error on 
each parameter is given by half the 16 — 84 percent percentile 
range of the likelihood distribution (this would be equivalent 
to the ±lcr range for a Gaussian distribution). 

We note that the stellar ages, as well as the other phys- 
ical parameters, are derived by fitting the galaxy spectra as 
observed and so refer to the galaxies at the time they are ob- 
served. Because of the bright and faint magnitude limits in 
our sample, there is a strong correlation between luminosity 
and redshift. If uncorrected, it may introduce systematic ef- 
fects in correlations between age and luminosity (or mass). 
To avoid this, we correct our measured ages so that they 
are relative to the present, rather than to the point of ob- 
servation by adding to the measured age for each galaxy 
the look-back time to the redshift at which it is observed. 
The metallicity is left unchanged. This look-back time varies 
from 0.07 to 2.64 Gyr over the redshift range covered by our 
sample. The corrections mainly affect the most luminous 
galaxies, which are found out to higher redshifts. All light- 
weighted ages quoted throughout this paper refer to z = 0. 

As described in Paper I, our constraints on metallic- 
ity and age are sensitive to the observational signal-to-noise 
ratio (S/N) of the spectra. A median S/N per pixel of at 
least 20 is required to constrain metallicity reliably. For this 
reason, we consider here only those galaxies with a me- 
dian S/N per pixel greater than 20. This cut reduces our 
sample to 26,003 high-concentration galaxies in the red- 
shift range 0.005 < z ^ 0.22. This is the same sample of 
'early-type' galaxies as analysed in section 3.3 of Paper I. 
As shown there, the S/N requirement biases the sample 
towards low-redshift, high-surface brightness galaxies, but 
it does not introduce any bias in the luminosity, colour, 
velocity-dispersion and index-strength distributions of the 
early- type galaxy sample. 

Our selection of early-type galaxies, based only on the 
concentration parameter, includes early-type spiral galax- 
ies, galaxies with emission lines and active galactic nu- 
clei (AGN). We can divide our sample i nto five subclasses 
accord ing to the classification given by iBrinchmann et all 
( 2004) on the basis of the emission-line properties of SDSS- 
DR2 galaxies. Following their notation, we label as (1) un- 
classifiable ('Un class.') those ga l axies which cannot be clas- 
sified using the iBaldwin et al.l <ll98ll) diagram, i.e. mostly 
galaxies with no or weak emission lines; (2) AGN those 
galaxies with a substantial AGN contribution to their 



2 This might not be entirely correct for USa and ITy^, which 
are suspected to depend on the a/F e ratio at high metallicity 
iThomas et al . 12004 Ikorn et all2005l) . However, we did not find 
any discrepancy between the metallicities and ages derived in- 
cluding or excluding HSa+H^a (see section 2.4.2 of Paper I). 



emission-line fluxes; (3) star forming ('SF') those galaxies 
with S/N> 3 in H/3, Ha, [OIII]5007, [NII]6564 emission lines 
and for which the contribution to Ha from AGN is less than 
1 percent; (4) composite ('C') those galaxies with S/N> 3 
in the same four lines, but for which the AGN contribu- 
tion to Ha luminosity can be up to 40 percent; (5) low-S/N 
star forming ('low-S/N SF') those galaxies with S/N> 2 in 
Ha but without AGN contribution to their spectra. The 
sample includes 10,982 unclassifiable galaxies, 7782 AGN 
and 3018 low-S/N SF galaxies. There are 2858 composite 
and 1362 star-forming galaxies, accounting for roughly 11 
and 5 percent of the total sample, respectively. As expected, 
most galaxies in the sample have specific star formation rates 
(SFR/M*) less than 10 -10 ' 6 yr _1 , which are in most cases 
consistent with zero. We note that the distribution in spe- 
cific star formation rates for the subclasses of C and SF 
galaxies shows a peak at 10~ 10 '° yr^ 1 . 

Fig. shows the distributions of several observational 
quantities of interest to us for this sample: the r-band ab- 
solute magnitude M r , the g — r colour, the logarithm of the 
velocity dispersion (in km s _1 ) and the Mg2 index strength. 
The magnitude M r and the colour g — r represent rest-frame 
quantities. Th ey are corrected for e volution following the 
prescription of (Bernardi et aD i2005t) . as described in Sec- 
tion !3.1.T1 Also the Mg2 index strength is corrected to z = 
assuming passive evolution, as described in Section [3.21 In 
each row, the solid histogram represents the distribution for 
each of the above subclasses (whose fractional contribution 
to the total sample is given in the leftmost panel), while 
the distribution for the sample as a whole is represented by 
the shaded grey histogram. The distributions for the un- 
classifiable, AGN and low-S/N SF galaxies agree well with 
the distributions for the sample as a whole. In contrast, C 
and especially SF galaxies tend to concentrate in the low- 
luminosity, blue, low-velocity-dispersion and low-Mg2 tails 
of the distributions. 

The stellar metallicities, ages and stellar masses derived 
in Paper I for the galaxies in the sample are shown in Fig. [5] 
(shaded grey histograms). The distribution in metallicity ex- 
tends from OAZq to 2 Zq with a peak around Z = 1.6 Zq. 
The distribution in r-band light-weighted age extends from 
2.5 to 12 Gyr with a peak around t r = 9 Gyr. The distri- 
bution in stellar mass extends from 10 10 to 10 12 Mq with 
a peak around M» = 10 11 Mq. The right-hand panels of 
Fig- HI show the distributions of the associated errors, com- 
puted as one half the 68 percent confidence ranges in the 
estimates of \og(Z/Z e ), log(t r /yr) and log (M»/M Q ). The 
typical uncertainty in both metallicity and age estimates is 
about ±0.1 dex, while stellar mass is constrained to better 
than ±0.1 dex for the majority of the galaxies in the sample. 
The same is true when considering Unclass., AGN and low- 
S/N SF galaxies only (solid histograms). The distributions 
for SF (dotted line) and C (dashed line) galaxies extend to 
lower stellar metallicities, ages and stellar masses than the 
bulk of the sample, with slightly higher errors in all three 
parameters. 

In addition to the effects of age, metallicity and stellar 
mass, we are also interested in the influence of the a/Fe ra- 
tio on the observed properties of early-type galaxies. This 
ratio can be empirically quantified by the rel ative stre ngths 
of M g- and Fe-based absorption-line indices jThomas et alJ 
2003). We use here Mgb/(Fe), where (Fe) is the average of 
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Figure 1. Distributions of the r-band absolute magnitude, g — r colour (^-corrected to z = and corrected for evolution following 
Bernardi et al. 2005), logarithm of velocity dispersion (in km s — 1 ) and Mg2 index strength (in mag) for a sample of 26,003 early-type 
(high-concentration) SDSS-DR2 galaxies with median S/N per pixel greater than 20 (grey histogra ms, repeated in all rows). F rom top to 
bottom, solid histograms in each row show the distributions for subsamples of galaxies classified by Brinchm ann et al.l J2004I) according 
to their emission-line properties as star-forming ('SF'), composite ('C'), low signal-to-noise star-forming ('low-S/N SF'), AGN ('AGN') 
and unclassifiable ('Unclass.', these are mainly galaxies with no emission lines). The fractional contributions by the different classes to 
the total sample of 26,003 early-type galaxies are indicated in the left-most panels. All histograms are normalised to unit area. 



the Fe5270 and Fe5335 index strengths. In what follows, we 
compare the observed Mgb/(Fe) ratio of a galaxy to that of 
the model that best reproduces the 5 spectral features men- 
tioned above. Since the models have scaled-solar abundance 
ratios, any discrepancy A(Mgb/(Fe)) between observed and 
model index strengths can then be interpreted as a depar- 
ture of the Q-elements-to-iron abundance ratio from solar in 
the observed galaxy. 

Working with the difference in index ratio rather than 
with the observed value of Mgb/(Fe) allows us to take 
into account the dependence of the index strength on 
age. It is also interesting to check whether the difference 
A(Mgb/(Fe)) traces linearly the true a/Fe abundance ra- 
tio over the whole parameter space covered by our galax- 
ies. It is not possible to check this self-consistently on our 
sample with our own models, since they do not include 
variations in element abundance ratios. For this test, we 
used Simple Stellar Population ( SSP) models w i th va riable 
element abundance ratios from iThomas et alJ i2003T) . We 
have considered SSPs with different a/Fe abundance ratios 
([a/Fe]=-0.2, 0, 0.2, 0.4, 0.6), with age varying from 1 to 13 
Gyr and metallicity varying from \og{Z / Zq) = —0.6 to 0.5. 
For each model we calculate the difference in Mgb/(Fe) with 
respect to the corresponding scaled-solar model at the same 
age and metallicity. We have looked at the relation between 
A(Mgb/ (Fe)) and the abundance ratio [a/Fe] at fixed metal- 
licity. We do not find significant differences in the slope and 



zero-point of the linear relations, averaged over age, for dif- 
ferent metallicities. Only at \og(Z / Zq) ~ —0.5 is the re- 
lation somewhat flatter; we note that this is not relevant 
for the galaxies in our sample, which predominantly have 
metallicities above \og{Z / Zq) = —0.2 (see Fig.|5J. The pro- 
portionality between A(Mgb/(Fe)) and [a/Fe] is confirmed 
by plotting A(Mgb/(Fe)) as a function of age and as a func- 
tion of metallicity for SSPs with different a/Fe ratio. At 
fixed [a/Fe], the measured A(Mgb/(Fe)) is constant over 
the age and metallicity ranges covered by our sample. These 
checks reassure us that the discrepancy A(Mgb/(Fe)) be- 
tween the observed and the best-fit model index strengths 
has the same proportionality with [a/Fe] over the age and 
metallicity ranges covered by the sample. 



3 PHYSICAL ORIGIN OF OBSERVED 

RELATIONS FOR EARLY- TYPE GALAXIES 

We use here the sample of 26,003 early-type galaxies de- 
scribed above to investigate the physical origin of two fun- 
damental relations obeyed by early-type galaxies: the colour- 
magnitude relation (Section 13. H and the Mg2-ov relation 
(Section 13.21 . We address the environmental dependence of 
these relations and of galaxy physical parameters in Sec- 
tion 13.31 We then explore in more detail how physical pa- 
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Figure 2. Left: distributions of the median-likelihood estimates 
of stellar metallicity (top), r-band light- weighted age (middle) 
and stellar mass (bottom) for the same sample of early-type galax- 
ies as in Fig. The dotted histogram shows the distribution for 
SF galaxies, the dashed histogram for C galaxies and the solid 
histogram for the rest of the sample (Unclass, AGN and low- 
S/N SF). Right: distributions of the associated errors, computed 
as one half the 68 percent confidence ranges in the estimates of 
log(Z/Z ), log(tr/yr) and log(M„/M Q ). 

rameters are related to the depth of a galaxy's potential well 
(Section 

3.1 The colour-magnitude relation 

The optical colours of early-type galaxies are strongly 
correlated with absolute magnitude, in the sense that 
bright g alaxies tend to be redder than faint galax- 
ies (e.g. iBauml Il959l: Ide Vaucouleursl Il96ll: iFaberl Il973t 
IVisvanathan fc Sandagdll977l: iBower et al.lll992l) . As men- 
tioned in Section Q the physical origin of this correlation is 
still a matter of controversy, mainly because of the lack of 
large samples of galaxies with accurate age and metallicity 
estimates. 

3.1.1 Observed colour-magnitude relation 

We consider here the colour-magnitude relation (CMR) de- 
fined by the rest-frame g — r colour and the absolute r-band 
magnitude M r . We compute these quantities using the SDSS 
model magnitudes, which are obtained from fits of the r- 
band surface brightness profile of each galaxy with either a 
de Vaucouleurs law or an exponential law (the best-fitting 
profile is also adopted in the other photometric bands). We 
A:-correct the observed g — r colour and the M r absolute mag- 
nitude to z — using the model in the spectral library of Sec- 
tion [5] that best reproduces the spectral absorption features 
of the galaxy, reddened to reproduce the observed colours at 



the redshift of the galaxy (this reddening correction is typ- 
ically small; see below). We also correct colours and mag- 
nitudes for evo l ution, adopting the correction estimated by 
iBernardi et alJ i2005l) from a sample of 39,320 SDSS early- 
type galaxies, i.e. we make the magnitudes fainter by 0.85z 
and the colours redder by 0.3z. From now on, we denote 
by M r and g — r the AH-e-corrected (but not dereddened) 
quantities. 

Fig. |3 shows the g — r, M r colour-magnitude relation 
defined in this way by early-type galaxies in the sample de- 
scribed in Section [5] We have arranged galaxies in bins of 
absolute magnitude and colour (0.15 and 0.01 width respec- 
tively). The grey scale indicates, for each magnitude bin, 
the relative distribution of galaxies in the different colour 
bins (we do not show bins containing less than 2 percent 
of the total number of galaxies at a given magnitude). The 
straight line in Fig. shows a 'robust' fit of these data, ob- 
tained by minimising the absolute deviation in colour as a 
function of magnitude. We note that this lies slightly to the 
blue of the ridge-line of the relation, reflecting the skewed 
colour distributions of Fig. Q The robust fit has a slope 3 of 
-0.024 ±0.002. This sl ope is consistent with that of -0.024 
obtained previously bv IBernardi et all <2005h 4 . The dashed 
lines in Fig.|3]show the mean positive and negative absolute 
deviations in colour (±0.050) relative to this fit. 

To better understand the influence of our galaxy selec- 
tion criteria on our results, we constructed a second sam- 
ple (we call this 'sample B') following the prescription by 
Bernardi et al.. We first select all objects classified as galax- 
ies, with measured velocity dispersions and with Petrosian r- 
band magnitudes (corrected for galactic reddening) between 
14.5 and 17.75 in the SDSS DR2. We then define as early- 
type those galaxies with f racDev_r > 0.8 and eclass < 0. 
We select in this way 59,907 early-type galaxies, for which 
we have g — r colours and M r absolute magnitudes. 5 

In Fig. 2] we compare the CMRs obtained from our pri- 
mary sample and from sa mple B in different r edshift bins 
(by analogy with fig. 2 of IBernardi et alJl2005l) . The black 
dot-dashed lines in Fig.|l]show robust fits to the relations for 
the two samples. The slope for sample B (—0.028 ± 0.001, 
upper panel) is in go od agreement with that obtained by 
IBernardi et alJ J2005L black solid line) and with that derived 
for our primary sample (—0.024 ± 0.002, lower panel). The 
colour-magnitude distributions of the galaxies in different 
redshift bins in sample B (up per panel) also agree reason- 
ably well with those in fig. 2 of Bernard^t^Ll J2005I) . Thus, 
our results agree with those of lB^rn^rdi^t^Ll 120051) under 
similar assumptions. 

We note that our primary sample differs from sam- 
ple B mainly through the inclusion of blue galaxies at 
intermediate-to-low luminosities. These are predominantly 
galaxies in the SF and C classes, which are likely to be 

3 The errors on the slopes of the fitted relations quoted through- 
out the paper have been estimated using a 'jackknife' method. 
The parameters of the fitted relations are summarised in Table HI 

4 We estimated this slope as (.CM^on I M i usi ng the values 
reported in tables 1 and 2 of IBernardi et all 820051 . see their ap- 
pendix Bl). ^^^^^^^^^^ 

5 Note that the sample used bv IBernardi et all i2005l) is slightly 
smaller than the SDSS DR2. Note also that sample B is larger 
than our primary sample, mainly because of our cut in S/N ratio. 
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bulge-dominated spiral galaxies with detectable emission 
lines (eclass>0, which excludes them from sample B). The 
contamination of our primary sample by these galaxies also 
results from the fact that we do not make any distinction 
on the basis of galaxy environment. While the slope of the 
CMR in low-density environments is not significantly dif- 
ferent from that found for cluster early-type galaxies, the 
scatter about the relation a ppears to increase from high- 
to low -density environments llLarson et alll 980 ; Hog g et alJ 
2004) . Also, at the m edian redshift z ~ 0.1 of our sample, 
IPimbblet et alJ J2002t) find that, at fixed luminosity, early- 
type galaxies tend to have bluer colours in t he outskirts than 
in th e centres of galaxy clusters (see also lAbraham et alJ 
Il99rf) . They interpret this as a difference in age of the galax- 
ies in different environments. 

As in other work on early- type galaxies (e.g. 
iBernardi et"afll2005l) . we do not include any correction for 
dust attenuation in the CMR presented here. We can esti- 
mate the influence of dust on the CMR, using the model 
spectral fits described in Section [5] The g- and r-band dust 
attenuation can be estimated from the difference between 
the emission- line corrected r — i fibre 6 colour of the galaxy 
and the r — i colour of the redshifted, dust-free model, as- 
suming a single power-law attenuation curve (see also Pa- 
per I). To quantify the effect of dust attenuation on the 
CMR, we show in Fig. |^| the average attenuation vectors 
(red arrows) for galaxies falling into different bins of g — r 
and M r (dashed boxes). As expected, the average attenu- 
ation is reduced if C and SF galaxies are excluded (blue 
arrows). In both cases, attenuation by dust becomes more 
important at lower luminosities. Correcting for dust atten- 
uation would thus steepen the relation. When de-reddening 
the colours and magnitudes of the galaxies in the sample ac- 
cording to the average attenuation vectors shown in Fig. 
we derive a slope of —0.035 ± 0.002. Tabic summarizes the 
parameters of the fits to the (fc-corrected) colour-magnitude 
relation, with and without evolution and dust corrections. 

3.1.2 Physical origin of the colour-magnitude relation 

Fig.|S|shows how metallicity, age, a/Fe ratio and stellar mass 
change along the colour-magnitude relation. We have binned 
the g — r, M r plane in the same way as in Fig. El the bright- 
ness of each g — r bin at a given M r being proportional to the 
number of galaxies falling into this bin (see Section 13.1.11 . 
In each panel, the colour code reflects the average properties 
of the galaxies falling into each colour-magnitude bin, as in- 
dicated. Panels (a)-(c) show how (the median-likelihood es- 
timates of) age, metallicity and stellar mass are distributed 
along the CMR. Panel (d) shows the distribution of the off- 
set A(Mgb/(Fe)) between observed and predicted Mgb/(Fe) 
index strengths, which traces the a-elements-to-iron abun- 
dance ratio (Section [2). 

The most remarkable result from Fig. [5] is that the 
colour-magnitude relation is primarily a sequence in the 
stellar mass of early-type galaxies (panel c; note that the 
gradient in stellar mass is not exactly parallel to the mag- 
nitude axis, reflecting a non-constant mass-to-light ratio in 

6 fibre magnitudes are obtained from the flux measured within 
an aperture of radius 1.5 arcsec, i.e. equal to the fibre radius. 
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Figure 3. Relation between (fc+e-corrected) g — r colour and r- 
band absolute magnitude for the same sample of early-type galax- 
ies as in Fig. The grey scale indicates, for each magnitude bin, 
the relative distribution of galaxies in the different colour bins 
(normalised along the colour axis). The solid line represents the 
relation between g — r and M r , fitted as described in the text. 
The dashed lines show the vertical scatter about the relation. 

the optical bands). Another striking result is that the aver- 
age metallicity increases from ~ 0.8Zq to ~ I.GZq from the 
faintest to the brightest galaxies along the relation, the aver- 
age light-weighted age increasing by less than 3 Gyrs (from 
6.5 and 9.0 Gyr) over the same interval of 5 magnitudes in 
M r (panels a and b). The marked increase in metallicity is 
accompanied by an increase in the a-elements-to-iron abun- 
dance ratio (panel d). This is consistent with the finding that 
giant elliptical gal axies have higher [Mg/Fe] than faint ellip- 
tical galaxies (e.g. | Wort hey et al"||l992|; [ Traeer et alJl2000al : 
iKuntschner et alJl200lHThomas et alJl2004l) . The gradients 
in both metallicity and age are not exactly parallel to the 
relation, but at fixed magnitude, bluer galaxies are more 
metal-poor and younger than redder galaxies. Thus, both 
metallicity and age contribute to the scatter about the rela- 
tion. 

It is of interest to compare the relative contributions by 
age and metallicity to the scatter about the CMR. We note 
that we expect dust attenuation and changes in the [a/Fe] 
ratio to contribute negligibly to the scatter, since the scatter 
in g — r colour of dust-free, [a/Fe]=0 models with the ages 
and metallicities of the galaxies in our sample (±0.041) 7 is 
very similar to the observed scatter (±0.050). For compari- 
son, the scatter in g — r colour of models with the metallic- 
ities of the galaxies in our sample but with a fixed imposed 
(average) age at fixed M r is only ±0.022, while that of mod- 
els with the ages of the galaxies in our sample but with a 
fixed imposed (average) metallicity at fixed M r is ±0.028. 
This indicates that age and metallicity contribute a similar 

7 This is estimated by adopting, for each galaxy in the sam- 
ple, the median of the likelihood distribution in g — r obtained 
by weighting each model in the library of Section [3] by \ 2 = 
[(log Z mod - logZ)/a logZ } 2 + [(logt mod - logt r )/cr log t r ] 2 . 
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Figure 6. Relation between (fc+e-corrected) g — r colour and M r absolute magnitude for the same sample of early-type galaxies as 
in Fig. Q The colour-magnitude relation has been binned and colour-coded to reflect the average stellar metallicity (panel a), r-band 
light- weighted age (panel b), stellar mass (panel c) and the o/Fe-indicator A(Mgb/(Fe)) (panel d) of the galaxies falling into each bin. 
The shading indicates the fraction of galaxies populating each colour bin at fixed magnitude, according to the map of Fig. El 



amount to the scatter in the colour-magnitude relation for 
early- type galaxies. 

Figure[7|better quantifies the intrinsic scatter in the two 
physical parameters and how much it determines the scatter 
in colour as a function of stellar mass. The left-hand panels 
show the distribution in light-weighted age (or metallicity) 
as a function of g — r colour in six bins of stellar mass (the 
median and the 16 — 84 interpercentile range are given by 
the red points and associated error bars). The right-hand 
panels show the marginalized distribution in age (metallic- 
ity) compared with a Gaussian of width equal to the aver- 
age error in age (metallicity) for the galaxies in each stellar 
mass bin (dotted curve). From the left-hand plot we can see 
that at masses above 10 11 Mq the scatter in light- weighted 
age can be entirely accounted for by the measurement er- 
rors, which are on average less than 0.1 dex. We can thus 
say that all massive ellipticals of given mass have the same 
mean age (within the errors). At lower masses the scatter 
in light- weighted age becomes larger than the typical error, 
and the distribution is skewed toward younger ages. At fixed 
mass, there is a correlation between age and colour, which 
saturates for the reddest galaxies. We also detect a small 
intrinsic scatter in the mass-metallicity relation in all mass 
bins (right-hand plot), which contributes to the scatter in 
the CMR at fixed mass, even at the high-mass end. 

It is important to check that aperture effects do not 
introduce any spurious trend in Fig. |S| This could arise 
because colour, luminosity and stellar mass pertain to the 



galaxy as a whole, while age, metallicity and A(Mgb/(Fe)) 
are measured within a fixed fibre aperture of radius 1.5 arc- 
sec, which samples different fractions of the total light in 
galaxies with different apparent sizes. In section 3.4 of Pa- 
per I, we already showed that aperture effects are not ex- 
pected to affect significantly the ages and metallicities de- 
rived for the galaxies in our sample. To probe the influ- 
ence of aperture effects on the results of Fig. |S| we have 
arranged galaxies in several bins of apparent size r$o,r, de- 
fined as the r-band Petrosian half-light radius. Galaxies with 



1.5 < 



< 2.4 arcsec are found out to redshift 0.2. 



We plotted the equivalent of Fig. |S] for each bin of rso,r 
within this range. In all cases, we found the same increase 
in metallicity, A(Mgb/(Fe)) and age along the relation as 
in Fig. EJ although the trends were noisier because of the 
smaller numbers of galaxies in individual r 5 o, r bins. In our 
sample, galaxies along the CMR have a roughly constant ap- 
parent size rso,r ~ 2.3 arcsec at all magnitudes fainter than 
M r = —23.5 (corresponding roughly to stellar masses less 
than lO n ' 5 M0). The brightest galaxies tend to have larger 
apparent sizes, but stellar metallicity, age and A(Mgb/(Fe)) 
do not show any significant trend with apparent size. Thus, 
we are confident that aperture effects are not responsible for 
the trends seen in Fig. 

For completeness, we repeated the analysis of Fig. fus- 
ing sample B in place of our primary sample (Section 13.1.11 . 
We found no significant difference in this case with respect 
to the results described above. We also repeated the analysis 
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Figure 7. Left panels: distribution in light-weighted age (left plot) and stellar metallicity (right plot) as a function of g — r colour for 
different bins of stellar mass. The red dots are the median age (metallicity) at fixed colour, while the error bars give the 68 percent 
percentile range in the age (metallicity) distribution. Right panels: histogram of light-weighted age (metallicity) compared to a Gaussian 
(dotted curve) of width equal to the average error on age (metallicity) in each stellar mass bin. In each panel the rms scatter and the 
average error in age (metallicity) are also given. 



using our primary sample, but correcting the observed CMR 
for the effects of dust attenuation (Section lbM.lH . Again, the 
main results in this case did not change. We further checked 
the influence of a possible contamination of the trends iden- 
tified in Fig. HO by late- type galaxies, in particular, with re- 
gard to the variation in age along the relation. Excluding SF 
and C galaxies (which tend to have younger ages than the 
bulk of the sample) had no significant effect on the results. 
In fact, C galaxies alone display the same trends as found in 
Fig. [S] for the sample as a whole. The average metallicity of 
SF galaxies also increases with luminosity along the CMR, 
although these tend to have light-weighted ages clustered 
around a mean value log(t r /yr) = 9.6 almost independently 
of luminosity. We note that SF galaxies represent only about 
5 percent of our primary sample, and they occupy g — r, M r 
bins where the fraction of galaxies at a given magnitude is 
less than 2 percent of the total sample (not displayed in 
Fig.©. 



In summary, we have shown that the colour-magnitude 
relation is primarily a sequence in galaxy stellar mass. Both 
the chemical composition (i.e. the total metallicity and the 
a-elements-to-iron abundance ratio) and the age of ellipti- 
cal galaxies depend mainly on stellar mass, increasing along 
the relation. At the high-mass end of the relation, the age 
spread is negligible and consistent with the errors. In this 
regime, the scatter in the CMR is determined by the small 
scatter in the mass-metallicity relation. At lower masses, the 
distribution in age becomes broader, with a spread toward 
younger ages, which correlates with colour and is thus the 
main contributor to the scatter about the colour-magnitude 
relation at the low luminosity end. 



3.2 The Mg 2 -0V relation 

In this section, we focus on another observational relation- 
ship between the stellar populations and the structural prop- 
erties of early-type galaxies: the relation betw een Mg2 index 
stren gth and galaxy velocity dispersion ay feender et alJ 
1993). The difficulty in interpreting this relation comes from 
the complex translation of Mg2 index strength and ay into 
physical parameters. 

The Mg2 index strengths considered here are corrected 
for broadening due to velocity dispersion. They are nor- 
malised to a fixed velocity dispersion of av = 200 km s~\ 
corresponding roughly to the average velocity dispersion 
of galaxies i n our sample. The norm alization is achieved 
by using the lBruzual fc Charlod Uo03) models to compute 
the difference in Mg2 index strength between the observed 
and reference velocity dispersion at the metallicity of each 
galaxy. 8 For a consistent comparison with the CMR and 
with the evolution-corrected light-weighted ages, we also 
correct the Mg2 index strengths so that they are relative 
to the present rather than to the redshift of observation, as- 
suming passive evolution. The corrections are obtained by 
fitting simple linear relations between Mg2 index strength 
and time for SSPs of different metallicities and velocity dis- 
persions. The Mg2 value measured for each galaxy is then 
corrected by adding the expected change in index strength 
over a time interval equal to the look-back time at the red- 
shift of the galaxy. 9 



8 This correction varies from —0.002 to 0.004 over the range in 
velocity dispersion covered by the sample. The average absolute 
correction is only 0.8 percent of the 5 — 95 percent percentile range 
in Mg2 index strength. 

9 The fitted slopes, averaged over velocity dispersion, are 
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Figure 4. Colour-magnitude relation for the same sample of 
early-type galaxies as in Fig. (lower panel) an d for a sample 
obtain ed using the same selection criteria as in Bcrnardi ct al. 
(2005) (upper panel). The colours and luminosities of the 
galaxies in both sam ples are corrected for evolution following 
IBernardi et al] <2005l i.e., making g — r redder by 0.3z and M r 
fainter by 0.852:). Different colours refer to different redshift in- 
tervals, 0.02 < z < 0.07 (cyan), 0.07 < z < 0.09 (blue), 
0.09 < z < 0.12 (green), 0.12 < z < 0.15 (yellow), 0.15 < z < 0.2 
(red). For each redshift bin, solid and dashed lines give the median 
and the 68 percent range in colour as a function of luminosity. 
Black dot-dashed lines show robust fits to the relations for the 
two samples (including galaxies at all redshifts), while the black 
sol id line in the upper p anel is the corresponding relation fitted 
by IBernardi et al] ]2005l) with a maximum-likelihood technique 
(see their fig. 2). 




1 



Figure 5. Relation between (fc+e-corrected) g — r colour and r- 
band absolute magnitude for the same sample of early-type galax- 
ies as in Fig. The red arrows indicate the average correction 
for dust attenuation on the colour and absolute magnitude of 
the galaxies falling into different g — r, M r bins (indicated by the 
dashed boxes) along the fitted relation. The blue arrows indicate 
the average dust correction when galaxies in the star-forming and 
composite classes are excluded. 
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Figure 8. Relation between Mg2 index strength (mag) and veloc- 
ity dispersion logov (km s _1 ) for the same sample of early- type 
galaxies as in Fig. As in Fig. |3| the grey scale indicates, for 
each velocity dispersion bin, the relative distribution of galaxies 
in the different Mg2 bins. 



Fig- IU shows the Mg2-o"v relation defined by the early- 
type galaxies in our primary sample. We have arranged 
galaxies in narrow bins of Mg2 and log cry (of widths 0.01 
and 0.036, respectively). By analogy with Fig. [3] the grey 
scale in Fig. |H| indicates, for each log cry- bin, the relative 
distribution of galaxies in the different Mg2 bins (we do not 
show bins containing less than 2 percent of the total number 
of galaxies at given velocity dispersion). The solid straight 
line shows the result of a robust linear fit obtained by min- 
imising the absolute deviation in Mg2 strength as a function 
of log ay. This has a slope of 0.25 ±0.01, i.e. very close to that 
found bv IBernardi et al ] J2003dl) for a sample of co-added, 
high-quality spectra of SPSS early-typ e galaxies (see also 
iGuzman et allll992l IColless et al.lll999l) . The dashed lines 
in Fig. |S] are the mean positive and negative absolute de- 
viations in Mg2 strength (±0.024 mag) relative to this fit 
(also in agreement with previous work). These results are 
summarized in Table Q We also indicate in the table the 
slightly lower slope of the Mg2-ov relation obtained for sam- 
ple B (0.19 ±0.002). This is because sample B includes fewer 
(emission-line) galaxies with weak Mg2 absorption than our 
primary sample at low ay and more (low-S/N) galaxies with 
weak Mg2 absorption at large ay (Section l3.1.11 . 

The tightness of the Mg2-ov relation has been used 
in the past to argue that early-type galaxies have nearly 
coeval stellar populations, perhaps within ~ 15 percent, 
at given co mposition (e. g. [Bend er et al.l ll993). However, 
as noted by iTraeer et all (l2000al) . the small scatter about 
the relation could also conceal significant age spreads, if 
these are accompanied by metallicity spreads such that the 
Mg2 strength rem ains roughly cons tant at fixed velocity dis- 
persion (see also Ijorgensenl fl999). Assumin g a moderate 
anticorrelation between age and metallicity. IColless et all 



0.0056,0.0044,0.0064,0.0077 mag 
-0.7, -0.4, 0, 0.4 respectively. 
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Figure 9. Distribution of residuals relative to the mean Mg2— cry 
relation fitted in Fig. |H] for different values of log (cry /km s —1 ) 
(indicated in each panel). The number of galaxies is normalised 
to the total in each bin of velocity dispersion. The dotted curve 
(repeated in each panel) is a Gaussian distribution of width equal 
to the average absolute deviation in the Mg2 -cry relation (see 
Table 0. 



(1999) find that the intrinsic scatter in the Mg2-cry rela- 
tion translates into an upper limit of 40 (50) percent on 
the age (metallicity) spread. Moreover, it has been pointed 
out that, at fixed cry, the distribution of Mg2 residuals 
relative to the mean Mg2-ov relation is not symmetric 



tends to increase at low velocity dispersion ( Burstcin et al 


1988; 


Bender et al.ll993l:lj0r£ 


rensen et al. 1996: Traeer et al 


2000a 


; IWorthev & CollobertJ 12003!) . Our data confirm this 



finding. Fig. [§] shows the distribution of Mg2 residuals at 
different velocity dispersions (increasing from the top-left to 
the bottom-right panels) for our primary sample. At large 
velocity dispersion, the distribution of residuals is symmet- 
ric and centred around zero, while at smaller velocity disper- 
sions, a tail of negative residuals appears. This could result 
from either an age or a metallicity spread, or both. The 
observation that galaxies with negative residuals are often 
morphologically disturbed (presumably because of recent 
star formation) led to the idea that low Mg2 values at fixed 
cry a re associated with younger ages llSchweizer fc Seitzerl 
Il992l) . 

The availability of independent constraints on the ages 
and metallicities of early-type galaxies in our sample allows 
us to re-examine the physical origin of the Mg2-cry relation. 
Fig. 1101 shows how metallicity, age, cv/Fe ratio and stellar 
mass change along this relation. We have binned and colour- 
coded the Mg2 , cry plane to reflect the average properties of 
the galaxies falling into each bin. As in Fig- [HI the brightness 
of each Mg2 bin at a given ay is proportional to the number 
of galaxies falling into this bin. Panels (a)-(c) show how (the 
median-likelihood estimates of) age, metallicity and stellar 
mass are distributed along the Mg2 -ay relation. Panel (d) 



shows the distribution of the offset A(Mgb/(Fe)) between 
observed and predicted Mgb/(Fe) index strengths, which 
traces the a/Fe ratio (Section 

The most remarkable result from Fig. 1101 is that, like 
the CMR, the Mg2-oy relation for early-type galaxies ap- 
pears to be primarily a sequence in stellar mass (panel c) . In 
fact, stellar mass correlates tightly with ve locity dispersion, 
which is a tracer of dynamical mass (e.g.. ICappellari et alJ 
2006). Panel (a) further shows that stellar metallicity in- 
creases along the relation, from ~ 0.8Zq to 1.6Zq from 
low to high velocity dispersions. Galaxies with large ay are 
also older than those with low ay, though for log cry < 2.3 
age appears to correlate with Mg2 index strength as well. 
Another striking result from Fig. 1101 is the similarity be- 
tween the metallicity and A(Mgb/(Fe)) gradients along the 
Mg2-ov relation (panels a and d). This shows that massive 
early-type galaxies with large velocity dispersions are both 
more metal-rich and more abundant in a elements relative 
to iron than less massive galaxies. 

As in the case of the CMR (Section 13.1.21 . we have 
checked that the trends identified above are not caused 
by aperture effects. In particular, we find similar trends in 
stellar mass, metallicity, age and A(Mgb/(Fe)) along the 
Mg2-ov relation when considering galaxies in narrow ranges 
of r"5o,r. Galaxies along the relation have roughly constant 
apparent size rso.r- ~ 2.3 arcsec, implying that the average 
fraction of total galaxy flux sampled by the fibre is almost 
constant. 

As above for the CMR, we can quantify the relative con- 
tribution of metallicity and age to the scatter in the Mg2-cry 
relation. We note that there must be another parameter, 
most likely the a-elements-to-iron abundance ratio, that is 
responsible for the scatter in the relation. In fact, the scatter 
in Mg2 of [a/Fe]=0 models with the ages and metallicities of 
the galaxies in our sample is only ±0.0118, i.e. much smaller 
than the observed scatter (0.0237). 10 We obtain a similar 
scatter (±0.0116) in the Mg2 index strength of models with 
the same metallicity as the galaxies in our sample but with 
a fixed imposed (average) age at fixed log ay. On the other 
hand the scatter in Mg2 index strength of models with the 
ages of the galaxies in our sample but with a fixed imposed 
(average) metallicity at fixed log ay is only ±0.007. This in- 
dicates that metallicity has a stronger influence than age on 
the scatter in the Mg2-ov relation. 

Figure 1111 is obtained in a similar way as Fig. |7| As in 
Fig. EI the left-hand panels of Fig. tTTI show that the scatter in 
light-weighted age of high-mass ellipticals (M* > lO n M ) 
is negligible and consistent with the errors on the age esti- 
mates. Only at masses below 1O 11 M0 does a significant tail 
of younger ages appear. These young galaxies are respon- 
sible for the scatter in Mg2 at fixed mass, but the relation 
between age and index strength saturates for Mg2 > 0.25. 
Also in agreement with Fig. the right-hand plot shows 
that we detect (above the measurement error) a small scat- 
ter in metallicity at fixed mass in all mass bins. There is a 
correlation between metallicity and Mg2 at fixed mass, in 



10 For each galaxy in the sample, we adopt here the median of the 
likelihood distribution in Mg2 obtained by weighting each model 
b yX 2 = [(log^m d-logZ)/cr logZ ] 2 + [(logt mod -logt r )/o- logtr ] 2 . 
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Figure 10. Relation between Mg2 index strength (mag) and velocity dispersion (km s — 1 ) for the same sample of early-type galaxies as 
in Fig. Q The relation has been binned and colour-coded to reflect the average stellar metallicity (panel a), r-band light- weighted age 
(panel b), stellar mass (panel c) and o/Fe-estimator A(Mgb/{Fc)) (panel d) of the galaxies falling into each bin. The shading indicates 
the fraction of galaxies populating each Mg2 bin a fixed log try- , according to the map of Fig.|U 




1.1 0.15 0.2 0.25 0.3 0.35 1 2 3 4 5 

^§3 frac. of galaxies/bin 



35 1 2 3 4 
frac. of galaxies/bin 



Figure 11. Left panels: light- weighted age (left plot) and stellar metallicity (right plot) as a function of Mg2 index strength (mag) in 
different bins of stellar mass. The red dots with error bars represent the median and the 68 percent percentile range of the age (metallicity) 
distribution at given Mg2 index strength. Right panels: histogram of light-weighted age (metallicity) compared to a Gaussian (dotted 
curve) of width given by the average error in age (metallicity) for the galaxies falling into each stellar mass bin. The scatter and the 
average error in age (metallicity) are also indicated. 
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Table 1. Parameters of the relations fitted by minimising the 
absolute vertical deviations. Sample A refers to our primary sam- 
ple, while sample B is obtained by selecting galaxies according to 
iBernardi et alj 120051) criteria. 



Sample 


Slope 


Intercept 


Scatter 


Colour-magnitude relation 


sample A a 


-0.015 ± 0.001 


0.40 


0.051 


sample B a 


-0.014 ±0.002 


0.41 


0.043 


sample A b 


-0.024 ± 0.002 


0.22 


0.050 


sample B 6 


-0.028 ± 0.001 


0.14 


0.043 


sample A c 


-0.025 ± 0.001 


0.14 


0.051 


sample B c 


-0.015 ± 0.002 


0.37 


0.043 


sample A d 


-0.035 ± 0.002 


-0.05 


0.050 


sample B d 


-0.029 ± 0.002 


0.09 


0.044 


Mg2— ov relation 


sample A e 


0.23 ±0.01 


-0.29 


0.024 


sample B e 


0.17 ±0.003 


-0.15 


0.023 


sample 


0.25 ±0.01 


-0.31 


0.024 


sample B^ 


0.19 ±0.00 


-0.18 


0.023 



"fc-corrected colour and magnitude. b k+ e-corrected colour 
an d magnitude (the evo lution correction is the one provided 
bv lBernardi et al"ll2005)) . c fc-corrected colour and magnitude 
dereddened applying the average dust corrections of Fig. [F] 
d Fully corrected colour and magnitude. e Mg2 index strength 
corrected for velocity dispersion. ^Mg2 index strength cor- 
rected for velocity dispersion and evolution. 

particular for galaxies with Mg2 > 0.25 (left-hand panels), 
where light-weighted age saturates to a constant mean value. 

The Mg2-av relation is, therefore, primarily a sequence 
in galaxy stellar mass. It reflects the fact that early-type 
galaxies form a sequence of increasing total stellar metal- 
licity and a-elements to iron abundance ratio from shal- 
low to deep potential wells. Our results also confirm the 
tren d of increasing age with increasing veloci t y dispersion 
(e.g. ICaldwell et al.ll20o3 iThomas etal]|2005l iNelan et all 
2005). At high masses, the small scatter in the relation cor- 
relates with the scatter in stellar metallicity at fixed stellar 
mass, while light-weighted age is almost independent of in- 
dex strength. At the low-mass end of the relation, the frac- 
tion of young ellipticals increases and variations in age be- 
tween ~ 4 and ~ 8 Gyr correlate with index strength at 
fixed mass. 



3.3 Environmental dependence 

Several authors have addressed the dependence on environ- 
ment of the properties of the stellar populations in early-type 
galaxies. Studies of the colour-magnitude relation and the 
relations between absorption indices and velocity dispersion 
on relatively small samples of early-type galaxies in different 
environments have shown that galaxies in low-density envi- 
ronments tend to be younger and m ore metal-rich than those 
in high-density envi r onments (e.g. |Kuirtschner_eJ_alJ [2002J; 
Tho mas et alJ 120051 : iDenicolo et alJ 12005!). I Bernardi et alJ 
( 1998), later confirmed bv lBernardi et al J ( 20061 on a larger 
sample of SDSS early-type galaxies, found differences in the 
Mg2— ov relation of galaxies in different environments, im- 
plying that galaxies in dense environments are at most 1 



Gyr older than galaxies in low-density environments and 
that they have the same metallicity. 

Kau ffmann et all i2004h provide estimates of environ- 
mental density for a sample of SDSS DR2 galaxies in the red- 
shift range 0.03 < z < 0.1 and with apparent r-band mag- 
nitude in the range 14.5 < r < 17.77, complete down to a 
stellar mass of 2 x 10 9 Mq. The density is expressed in terms 
of the number of spectroscopically-observed neighbouring 
galaxies (down to a fixed absolute magnitude) within 2 Mpc 
of projected radius and ±500 km/s in velocity difference 
from the target galaxy, corrected for galaxies missed due to 
fibre collisions (N ne i g h). We take advantage of these density 
estimates to address any environmental dependence of the 
physical properties of the stellar populations for the galaxy 
sample studied here. This can be achieved for only 1765 
galaxies in our sample, for which an estimate of N nc i g h is 
available. We consider three bins in environmental density, 
defined by Nneigh < 4, 4 < N ncig h < 7 and N ne igh > 7, 
which contain, respectively, 693, 388 and 684 galaxies. As 
mentioned in Section [5] we can classify the galaxies in our 
sample on the basis of their emission-line properties. As ex- 
pected, the sample is dominated by 'unclassifiable' galax- 
ies (without emission lines), but there is also contamination 
by galaxies with a low level of star formation (SF and C 
galaxies). Fig. 1121 illustrates the fraction of unclassifiable, 
SF, low-S/N SF, C and AGN galaxies as a function of envi- 
ronment. This plot quantifies the statement of Section l3.1.1l 
that the fraction of galaxies showing emission lines in our 
sample increases in lower-density environments. This class 
of galaxies also contributes to increase the scatter blueward 
of the colour-magnitude relation. 



In Fig. 1131 we explore how the CMR (left-hand panels) 
and the Mg2— ov relation (right-hand panels) depend on en- 
vironment. The relations found in the highest-density bin are 
compared to those defined in the lowest-density bin. The re- 
sults of the linear fits (also for the intermediate bin of N ne i g h) 
are given in Tabled Fig. 1131 and Tabled show that there is 
no systematic variation in the slope of the CMR as a func- 
tion of environment, while the Mg2— ov relation appears to 
steepen at low densities because of a larger fraction of galax- 
ies with low Mg2 index strength at low velocity dispersions 
(columns 3 and 6). Between the two extreme density bins 
there are differences of 0.006 ± 0.003 and 0.007 ± 0.003 in 
the zero-points of the CMR and the Mg2— ov relation, re- 
spectively (columns 4 and 7). This is in agre ement with the 
small shift of 0.007± 0.002 mag measured bv lBernardi et al] 
(1998) in the Mg2— ov relation of a sample of 931 early- 
type galaxies in different environments. We also identify a 
systematic increase of the scatter about both relations from 
high- to low-density environments, in agreement wit h earlier 
findings, as mentioned above (e.g. lHogg et al]l2004t) . 



It is of interest to understand how the change in the 
scatter along the two scaling relations reflects differences 
in the physical parameters of galaxies in different environ- 
ments. We note that the distribution in stellar mass does not 
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vary significantly 11 with environment, but that the median 
stellar mass of galaxies in the lowest-density bin is lower by 
about 0.05 dex than the median stellar mass in the highest- 
density bin (it increases from 8 x f0 10 to 9 x 10 10 M o ). Since 
stellar metallicity, age and element abundance ratio all in- 
crease with stellar mass, any effect induced by changes in 
stellar mass must be removed when quantifying variations 
in these parameters with environment. To do this, we cal- 
culate the median stellar metallicity, light-weighted age and 
element abundance ratio as a function of stellar mass for the 
sample as a whole (see Fig. I17H . For the 1765 galaxies with 
an estimate of environmental density, we then consider the 
offsets in \og(Z/Z Q ), log(t r /{/r), A(Mgb/(Fe)) from the me- 
dian values of these parameters at fixed stellar mass in the 
whole sample. The distributions in A[logZ] and A[logt r ] 
are skewed toward negative values, independently of envi- 
ronment, and this effect is stronger at small masses. If galax- 
ies at low densities are distributed preferentially to smaller 
masses than galaxies in dense environments, the distribu- 
tions in A[logZ] and A[logt r ] will show a stronger tail to- 
ward negative values in the low-density bin. To separate this 
effect from an intrinsic dependence of metallicity and age on 
environment, we further distinguish between galaxies with 
stellar masses above and below 1O 11 M0. 

The result of this analysis is shown in Fig. 1141 for the 
same three environmental bins as considered above. The 
distributions of the offsets in \og(Z / Zq), log(t r /yr) and 
A(Mgb/(Fe}) for the two low-density bins (N no i gn < 4 and 
4 ^ N nc igh < 7, solid lines) are compared to the correspond- 
ing ones in the high-density bin (N no i g h 7, dotted line in 
each panel), for massive and low-mass galaxies separately 
(red and blue histograms, respectively). The comparison 
is also summarized in Table [3] where we give the differ- 
ence of the average parameter offset between the highest- 
and the lowest-density bins, for galaxies with stellar mass 
above and below 10 Mq separately. There, we also indicate 
the probability for the two distributions to be drawn from 
the same parent distribution, according to a Kolmogorov- 
Smirnov test. Element abundance ratio, as expressed by 
A(Mgb/(Fe)), does not show any significant variation with 
environment, either in the average value or in the scatter. 
In contrast, light-weighted age (independently of mass) and 
metallicity (for massive galaxies) show a small dependence 
on environment in the sense that there is a higher fraction 
of young, metal-poor galaxies in low-density environments. 
Massive early-type galaxies in dense environments tend to 
be ~ 0.03 dex more metal-rich than their field counterparts. 
Similarly, the light-weighted ages of galaxies in dense en- 
vironments are ~ 0.02 dex older than in the field. It is 
interesting to mention that we also find a systematic in- 
crease in the scatter of metallicity and age from dense to 
low-density environments. From the highest- to the lowest- 
density bins the scatter in both metallicity and age increases 
by about 0.02 dex for massive galaxies (from 0.118 to 0.135 
for log(Z/Z Q ), from 0.083 to 0.103 for log(t r /yr)) and by 



1 When comparing the stellar mass distribution of galaxies in 
the low-density bin with that of galaxies in the high-density bin, 
the probability obtained from a Kolmogorov-Smirnov test is not 
low enough to reject the hypothesis that the two distributions are 
drawn from the same parent distribution. 




Figure 12. Fraction of unclassifiable ('Unclass'), star-forming 
('SF'), low-S/N star-forming ('low-S/N SF'), composite ('C') and 
AGN galaxies as a function of environment. Three different bins of 
environmental density are considered here, as given by N ne i gn in 
the plot. Each histogram is normalised to the number of galaxies 
in the corresponding environmental bin. While the high-density 
environments are strongly dominated by galaxies without emis- 
sion lines (unclassifiable), the fraction of star-forming, composite 
and AGN galaxies increases in lower-density environments. 



Table 3. Difference between the average parameter offset be- 
tween the N^gh < 4 and N ne ; gn 7 environmental bins, distin- 
guishing between galaxies with stellar masses above and below 
1O 11 M0. For each parameter, the second row gives the probabil- 
ity from a Kolmogorov-Smirnov test that the two distributions 
are consistent with being drawn from the same parent distribu- 
tion (we consider differences to be statistically significant only if 
the KS probability is below 1 percent). 





M, > ICFMq 


M, < 10 n M Q 


A[logZ] 


0.032 ± 0.006 


0.029 ± 0.006 




0.0030 


0.0217 


A[logt r ] 


0.019 ±0.005 


0.023 ± 0.006 




0.0053 


0.0017 


A[A(Mgb/(Fe»] 


0.008 ± 0.009 


0.024 ± 0.009 




0.5321 


0.1499 



about 0.01 dex for M* < lO n M (from 0.16 to 0.17 for 
log(Z/Z @ ) and from 0.12 to 0.136 for \og(t r /yr)). Although 
very small, these trends hint at a possibly very relevant en- 
vironmental dependence of metallicity and age. Future anal- 
ysis of larger samples of galaxies (provided, e.g., by the com- 
plete SDSS) with well characterized environmental proper- 
ties will allow to draw firmer conclusions. 

3.4 Correlations between physical parameters 

So far, we have considered separately the constraints set 
by the colour-magnitude and the Mg2-ov relations on the 
ages, chemical compositions and stellar masses of early-type 
galaxies in our sample. Here, we explore in more detail the 
potential correlations between these different physical pa- 
rameters, the dynamical mass and the stellar surface mass 
density of galaxies. 
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Figure 13. Colour-magnitude (left plot) and Mg2— cry relations (right plot) as a function of environment. Galaxies are divided into two 
disjoint bins of environmental density (increasing from bottom to top as indicated by N^gh in each panel). The linear relation fitted 
for each subsample and the scatter about it are shown in each panel by the black solid and dashed lines. The grey lines reproduce the 
relation fitted in the highest density bin. 



Table 2. Parameters of the colour-magnitude relation (columns 3,4,5) and Mg2— cry relation (columns 6,7,8) for a subsample of 1765 
galaxies with environmental estimates. The relations are fitted separately for three different bins of environmental density (given by 
Nnciglu column 1). The number of galaxies in each bin is given in column 2. The intercepts are given at M r = —21.5 and log cry = 2.25, 
which correspond roughly to the average magnitude and velocity dispersion of this subsample of galaxies. 







colour- magnitude 






Mg2 - T V 




Environment 


N gal 


Slope 


Intercept 


Scatter 


Slope 


Intercept 


Scatter 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Nncigh <C 4 

4 < N ncigh < 7 

Nncigh ^ 7 


693 
388 
684 


-0.021 ±0.004 
-0.026 ± 0.006 
-0.022 ± 0.003 


0.755 ± 0.002 
0.757 ± 0.004 
0.761 ± 0.003 


0.043 
0.041 
0.036 


0.27 ±0.03 
0.26 ±0.01 
0.24 ±0.02 


0.244 ±0.002 
0.242 ±0.001 
0.251 ±0.003 


0.023 
0.020 
0.019 



We first turn our attention to the relation between age, 
luminosity, velocity dispersion and stellar mass inferred from 
the distribution in the full parameter space described by the 
colour-magnitude and the Mg2-uv relations. Figs [(Jd and 
IlUb taken together indicate that brighter galaxies or galax- 
ies with large velocity dispersions are on average older than 
fainter or smaller velocity dispersion ones (c onsistent with 
the 'downsizing' scenario; ICowie et al.lll996l . see also Sec- 
tion^. In Fig. 115b . we show the relation between M r ab- 
solute magnitude, velocity dispersion cry and light-weighted 
age for the galaxies in our primary sample. At fixed luminos- 
ity, galaxies with large cry tend to be ol der than those with 
small cry (a result already poin ted out bv lForbes fc Ponmanl 
119991 and lBernardi et aljEooil) . Conversely, at fixed veloc- 
ity dispersion, luminous galaxies tend to be younger than 
faint galaxies. The dispersion in the M r -ov relation and 
the dependence of luminosity on age at fixed cry cannot be 
accounted for entirely by the fact that the stellar mass-to- 
light ratio M*/L r of an evolving stellar population increases 
with light- weighted age. This is shown in Fig. 115b . where we 
plot the relation between stellar mass and velocity disper- 
sion for the same galaxies as in Fig. 115b . The scatter of 



0.252 dex about this relation is consistent with the scatter 
of 0.579 mag in the relation between M r and cry. There ap- 
pears to be a real dependence of age on stellar mass at fixed 
velocity dispersion, in the sense that galaxies with large stel- 
lar masses tend to be slightly younger than those with low 
stellar masses (Table 2] provides simple linear fits to the re- 
lations shown in Figs 115b and b). 

The results of Fig. 1151 motivate us to investigate more 
closely the relation between stellar mass and dynamical mass 
for early- type galaxies. We can estimate the dynamical mass 
(including stars and dark matter) within the optical radius 
of a galaxy f rom a combination of velocity dispersion and 
optical size dCappellari et alJl200d) . Here we use the ra- 
dius containing 50 percent of the Petrosian flux in the r- 
band (Rso, r ) rather than the effective de Vaucouleurs radius 
(rdev), because we do not explicitly select early- type galax- 
ies on the basis of the shape of their light profile. Thus, a de 
Vaucouleurs profile may not always be an optimal fit to the 
photometric data, and the empirical half-light radius Rso,r is 
more straightforward to interpret. 12 To estimate the dynam- 

12 The de Vaucouleurs effective radii of the galaxies in our sam- 
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Table 4. Linear fits to the relations plotted in Figs ll5l and ll7l 





Slope 


Intercept 


Scatter 


log ay vs M r 
log ay vs log M* 


-0.125 ±0.004 
0.286 ± 0.020 


-0.464 
-0.895 


0.073 
0.071 


log t r vs log Mdyn 
log t r vs log M* 


0.115 ±0.056 
0.112 ±0.098 


8.628 
8.678 


0.085 
0.088 


log Z VS log M d yn 

log Z vs log M* 


0.164 ±0.010 
0.168 ±0.013 


-1.731 
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0.119 
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A(Mgb/(Fe)) vs \ogM dyn 
A(Mgb/(Fe)) vs logM, 
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Figure 14. Stellar metallicity, age and element abundance ra- 
tio (expressed as A(Mgb/(Fe))) as a function of environment 
for the subsample of 1765 galaxies with environmental esti- 
mates. We have estimated the offset in \og{Z I Zq), log(t r /yr) 
and A(Mgb/{Fe)) from the median value at fixed stellar mass 
calculated for the sample as a whole. The distributions in the off- 
sets are more skewed toward negative values at smaller masses. 
Thus, to further remove any effect induced by an environmen- 
tal dependence of the mass, we analysed separately galaxies with 
masses above and below lO n M0 (shown in red and blue respec- 
tively). For each stellar mass bin, the distributions in the off- 
sets are shown for galaxies with different environmental density 
(N ne igh < 4 in the lower panels and 4 ^ N ncigll < 7 in the up- 
per panels, solid line). Each distribution is compared to that for 
galaxies with Nneigh ^ 7, shown by the dotted line in each panel. 




-19 -20 -21 -22 -23 -24 9.0 9.5 1 0.0 1 0.5 1 1 .0 1 1 .5 1 2.C 

M -5log 10 h 70 log (M,/V G ) 



Figure 15. Relation between velocity dispersion and r-band ab- 
solute magnitude (panel a), colour-coded to reflect the average 
light- weighted age of the galaxies falling into each log ay— M T bin. 
For M r brighter than —20, lines of constant age are approximately 
parallel to the relation. Panel (b) shows the result of substituting 
absolute magnitude with stellar mass. 



ical mass within Rso,r, we must also account for the fact that 
the velocity dispersion of SDSS galaxies i s measured withi n 
a fixed 1.5 arcsec fibre radius. Following I.Torgensenl |l999), 
we assume that the radial profile of the velocity dispersion 
has a slope —0.04. The velocity dispersion at the radius Rso, r 
is then given by 



0"v(R5O,j 



' fibre 
R50.r 



and the dynamical mass within Rso,r by 



dy ri 



R50,i 



Mr. 



G 



(1) 



(2) 



pie are on average ~ 1.4 times larger than the Petrosian half-light 
radii. This arises presumably from the fact that the Petrosian flux 
represents typic ally 80 percent of the tota l flux in the de Vau- 
couleurs model iPadmana bhan et alj|2004l) . but this could also 
result (at least for a class of objects) from forcing a de Vau- 
couleurs fit to a flatter galaxy light profile. 
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Table 5. Correlation between stellar and dynamical mass. The 
first two columns give the radius within which the dynamical 
mass is estimated and the slope of the velocity dispersion profile 
assumed to correct for aperture effects. 



Figure 16. Relation between stellar mass and dynamical mass 
M-dyn estimated within the r-band Petrosian half-light radius (see 
equation |5J . Lines of different colours represent the median re- 
lations in different bins of light-weighted age, increasing from 
log(t, /yr) = 9.5 (purple) to 10.1 (yellow). The ratio between 
stellar and dynamical mass decreases with mass, as highlighted 
by the comparison with the one-to-one relation (dashed line). 



In this expression, G is the gravitational constant, and £ is a 
scaling factor given in Table|^| which depends on the shape 
of the velocity dispersion profile (this factor would be equal 
to 5 if we estimated the dy namical mass within th e effective 
de Vaucouleurs radius; see ICappellari et al]|200fj) . 

In Fig. 1161 we compare the dynamical masses obtained 
in this way to the stellar masses derived in Paper I for the 
early- type galaxies in our primary sample. 13 We focus on the 
slope of the relation, because the zero point is sensitive to 
our assumptions about the radius within which the dynami- 
cal mass is estimated, the conversion factor between dynam- 
ical and virial masses, and the stellar initial mass function 
(IMF; which we assume constant for all the galaxies). A lin- 
ear fit to the relation in Fig. 1161 yields 



(M./M ) oc (M d , n /M )(°- 783±0 - 019 ) 



(3) 



with a scatter of 0.13 dex. The specific correction applied 
to o v has little effect on the fitted slope. As a check, we 
tried steeper velocity dispersion pr ofiles, with slopes —0.0 6 
jMehlert et alJl2003f) and -0.066 llCappellari et al.ll2006l) . 
These yielded relations between log M» and log Md yn with 
slopes 0.785 ±0.021 and 0.786 ±0.017, respectively, i.e. con- 
sistent within lcr with the slope of the relation in Fig. 1161 
We also tested the effect of estimating the dynamical mass 
within the effective de Vaucouleurs radius rdev instead of 
R,50, r - The dynamical masses obtained in this case were sys- 
tematically higher, by ~0.13 dex, than those derived within 
R,50, r - The relation between log M* and log Md yn had a slope 



13 We recall that the total stellar mass is obtained by multiply- 
ing the stellar mass-to-light ratio derived from a fit of the fibre 
spectrum by the total luminosity of the galaxy. This assumes that 
the stellar mass-to-light ratio outside the fibre is the same as that 
inside the fibre (see Paper I). 
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Slope 


Intercept 


Scatter 


R-50,r 


-0.040 


7.14 


0.783 ± 0.019 


2.19 


0.127 


R-50,r 


-0.060 


6.98 


0.785 ±0.021 


2.18 


0.128 


R-50,r 


-0.066 


7.02 


0.786 ±0.017 


2.17 


0.128 


rdeV 


-0.040 


5.00 


0.808 ± 0.026 


1.93 


0.113 



of 0.808 ±0.026, i.e. slightly higher but consistent within la 
with the slope of the relation in Fig. 1161 These results are 
summarised in Table 

A robust result from Fig. 1161 therefore, is that the ra- 
tio between dynamical mass and stellar mass increases from 
the least massive to the most massive early-type galaxies in 
our sample. Structural non-homology does not appear to be 
responsible for this effect. We have built different subsam- 
ples of galaxies, based on the value of the Sersic index n 
fitted to the light profile in the SDSS database. The slope 
of the relation between log M* and log Mdyn within each 
subsample remains close to that of the relation in Fig. 1161 
changing from 0.847 for n — 3 to 0.801 for n = 5.5 (a de 
Vaucouleurs profile corresponding to n = 4). Instead, the 
decrease in M*/Md yn with stellar mass in Fig. 1161 is con- 
sistent with the increase in the dynamical mass-to-light ra- 
tio (Mdyn/L) of early-type galaxies implied by the Funda- 
mental Plane under the assumption of s tructural homology 
jBender et alJll992t iPierini et alJliool IZibetti et alJl200a 
see also ICaOTeUarFet^i . 2006, where no assumption on ho- 
mology is made). This is also consiste nt with the increase 
in Md y n/L with luminosity found by Padmanabhan et alJ 
( 2004) for a sample of 29,469 SDSS elliptical galaxies. The 
decrease of M*/Md y n with stellar mass could result from 
a more efficient mixing of dark matter and stars within 
the optical radius of massive galaxies relative to low-mass 
galaxies, as expected if the most massive early-type galax- 
ies assembled through multi ple mergers of dissipationless 
systems (see for discussio n IWjiitej |l980j; |Hmnj3hj^L£t^L 
2005HDe Lucia et al.l200rj:lBovlan-Kolchin. Ma fc Quataert 



20051) . 

The trend in M„ /Md yn with stellar mass shows a weak 
dependence on galaxy light- weighted age. This is shown in 
Fig. 1161 where lines of different colours indicate the median 
stellar mass as a function of dynamical mass for galaxies in 
various age bins, from log(i r /yr) = 9.5 (purple) to 10.1 (yel- 
low). Lines of constant age run parallel to the relation and, 
in spite of the small scatter, it appears that, at given dynam- 
ical mass, galaxies with more mass in stars are younger than 
those with small stellar mass (see also Fig. 115b ). This weak 
trend cannot be accounted for entirely by the larger amount 
of mass returned to the interstellar medium by evo lved stars 
in old er galaxies relative to younger ones. For the IChabrierl 
( 2003) IMF adopted here, the returned stellar mass fraction 
of a simple stellar population increases by about 0.03 dex 
from log(t r /yr) = 9.5 to 10 (with little depend ence on metal- 
licity; the differential change is similar for a lSalpeterfli"955l 
IMF). This effect can thus account for only about 10 per- 
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cent of the trend in Mt,/M dyn with age in Fig. ED 14 The 
bulk of the trend might result from a systematically higher 
baryonic fraction and/or higher efficiency of conversion of 
baryons into stars in young early-type galaxies relative to 
old ones. For example, if many of our early-types form by a 
merger of star-forming systems, those which currently have 
the youngest populations are presumably the most recently 
merged and so spent the longest time in the star-forming 
phase. 

We now examine in more detail how age, stellar metal- 
licity and a/Fe ratio depend on stellar and dynamical mass. 
This is shown in Fig. 1171 for the early-type galaxies in our 
primary sample (Table |3] provides simple linear fits to the 
relations shown in the figure). The relations followed by 
age, metallicity and a/Fe ratio as a function of stellar mass 
(right-hand panels) reflect the conclusions drawn from our 
analysis of the CMR and Mg2-ov relations in Section 13.11 
and 13.21 In particular, we find that light- weighted age in- 
creases (with a small scatter) with stellar mass in galaxies 
more massive than 1O 11 M0, while there is a clear indica- 
tion of a tail towards younger ages in less massive galax- 
ies (Panel d; see also fig. 12 of Paper I). This confirms the 
results of several previous studies of smaller samples of 



Caldwell & Rose 1998: Poggianti et al.1 2001: 


Thomas et alJ 


2005 


Ivan Dokkum & Ellis! 120031: iTreu et alJ 


2005). Stellar 



metallicity increases all the way from the least massive to 
the most massive galaxies in our sample (Panel e). The re- 
lation tends to steepen at stellar masses less than about 
~ 3 x 1O 1O M0, i.e. the characteristic mass scale of sev- 
eral observed bimodaliti es (see Paper I, iKauffmann et ail 
120031: iBaldrv et al-lEool) . We note that the dependence of 
stellar metallicity on stellar mass found here for early-type 
galaxies resemb l es qu ite closely the relation presented by 
iTremonti et alJ i2004T) for a sample of SDSS star-forming 
galaxies. 

The increase in metallicity with stellar mass, and hence 
the CMR, is a natural prediction of different scenarios of 
early-type galaxy formation. In models in which the galax- 
ies form monolithically in a single dissipationless collapse, 
star formation and chemical enrichment are generally as- 
sumed to be interrupted by the onset of galactic winds fol- 
lowing the major episode of star formation. Since massive 
galaxies with deep potential wells are able to retain their 
gas for much longer times and so to reach higher metallici- 
ties than less massive galaxies, t his scenario accounts for the 
CMR as a metallicity sequ ence llLarsorj[l97^; Bressan et alJ 



19941: iTantalo et alTll996l: lArimoto fe Yoshiilll987h . Early 



models of hierarchical galaxy formation neglected chemi- 
cal enrichment and faile d to produce sufficiently red, lumi- 
nous elliptical galaxies (IKauffmann et all Il993t ICole et all 
Il994f) . However, subsequent renditions of these models have 
been able to reproduce the observed slope and scatter of 
the CMR and of the Mg2-av relation by including chem- 
ical evolution and strong feedback even in massive galax- 
ies toget her with metallicity-dependent population synthesi s 
models (IKauffmann fc Charlotlll998t foe Lucia et alJl2004l) . 



14 For example, at M dyn ~ 1O 11 M0, stellar mass increases from 
log M t /M Q ~ 10.75 for log(i r /yr) = 10 to ~ 10.9 for log(t r /yr) = 
9.5. 



In this respect, the CMR of early-type galaxies does not 
represent a powerful tool to discrimina te between the mono - 
lithic and merger formation scenarios (iKavirai et al.l f2005'l. 

Additional clues about the formation of early-type 
galaxies may lie in their chemical abundance patterns. 
Fig. II 71 : shows that the a/Fe abundance ratio, as traced 
by the quantity A(Mgb/(Fe)) (Section[2J|, increases linearly 
with logM*. This result, which was anticipated in several 
previous studies of early-type galaxies (: Worthev et alJll992l : 
Ij0rgensenlll999l : lGreggidll997t) . has been quantified using 
stell ar population models with variable metal abundance ra- 
tios |Trager et alJl2000bl:lMehlert et alJl2003l:lThomas et alJ 
120051) . If we consider as a-enhanced those galaxies with 
A(Mgb/(Fe)) > 0.2 (to account for the typical observa- 
tional error on Mgb/(Fe)), Fig. I17t ' indicates that early- 
type galaxies with solar abundance ratios are found only 
at M* < 10 10 ' 5 Mq, corresponding roughly to velocity dis- 
persions ay < 100 km s -1 . This is similar to the conclusion 
drawn bv lKuntschner et alJ <l200ll) . based on the analysis of 
a sample of 72 early-type galaxies in groups and clusters. 

The observed range in a/Fe abundance ratio and the in- 
crease of this ratio with stellar mass constitute a challenge 
for galaxy formation models. If the enrichment in a elements 
relative to iron reflects the star formation timescale, as is 
assumed in standard chemical evolution models, the values 
A(Mgb/(Fe)) > 0.2 found in high-mass early-type galax- 
ies require star formation timescales of the order 1-2 Gyr 
(see Section 0J. Such short timescales are plausible for the 
onset of galactic winds in the classical monolithic collapse 
scenario. Hierarchical models can also produce star forma- 
tion histories which are peaked at high redshift for mas- 
sive early-type galaxies, although subsequent st ar formation 
is expected at lower redshift in most models dKauffmannl 
119961: iThomad 119991) . A possible solution to this problem 
may be the suppression of late gas cooling (and hence star 
formation) by AGN-driven outflows in massive early-type 
galaxies. As shown by various recent models of hierarchical 
galaxy formation accounting for the combined effects of star 
formation and black hole accretion, the feedback provided 
by active g alactic nuclei can halt star formation on short 
timescales (iGranato et al.ll200ll 120041: ISpringel et alJl2005l : 
ICroton et al.ll2006l : iDe Lucia et alJl200d) . This mechanism 
is expected to be most effective in massive halos hosting 
supermassive black holes, for which the timescale to driv e 
outflows could be as short as 1 Gyr (iGranato et al.ll2004h . 
Hence, the increase in the a/Fe abundance ratio with stel- 
lar mass adds important constraints to models of early-type 
galaxy formation. 

Fig. 1171 also allows us to compare how age, metallic- 
ity and a/Fe ratio vary with dynamical mass (left-hand 
panels) versus stellar mass (right-hand panels). In the left- 
hand panels, lines of different colours show the median 
relations between each parameter and Md yn followed by 
galaxies in different bins of stellar mass, increasing from 
log(M„/Af Q ) = 9.6 to 11.7 in steps of ~0.3 dex (from blue to 
red). In the right-hand panels, analogous lines show the me- 
dian relations between each parameter and M» followed by 
galaxies in different bins of dynamical mass. Fig. ll7H shows 
that, at fixed dynamical mass, galaxies with large M* tend 
to be younger than those with small M* (in agreement with 
the result of Fig. 115b above). This trend is noticeable for 
dynamical masses in the range 10 10 < M dyn < 1O U 5 M . 



18 A. Gallazzi, S. Chariot, J. Brinchmann, S.D.M. White 




9.5 10.0 10.5 11.0 11.5 12.0 9.5 10.0 10.5 11.0 11.5 12.0 

log(M dy „/M a ) log(M./M s ) 



Figure 17. Left: Correlations between different physical parameters and the dynamical mass estimated within the r-band Pet- 
rosian half-light radius (see equation [5J: (a) Light-weighted age; (b) stellar metallicity; (c) o/Fe-estimator A(Mgb/(Fe)). The 
grey scale indicates the density of points at each location in the diagram, normalised to the total number of galaxies at fixed 
log(Md an /MQ). The solid lines represent the median relations for galaxies in different bins of stellar mass, centred on log(M*/MQ) = 
9.6, 9.9, 10.2, fO. 5, 10.8, 11.1, 11.4, 11.7 (increasing from blue to red). Right: Correlations between the same physical parameters and 
stellar mass. In this case, the solid lines represent the median relations for galaxies in different bins of dynamical mass, centred on 
\og(M dyn /M Q ) = 9.6, 9.9, 10.2, 10.5, 10.8, 11.1, 11.4, 11.7, 12 (increasing from blue to red). 



In contrast, at fixed dynamical mass, stellar metallicity is 
almost independent of M t (Fig. 117^ ). In Fig. I17L galaxies 
with large M» tend to have slightly lower A(Mgb/(Fe)) and 
hence a/Fe than those with small M* at fixed Md yn . This 
is consistent with the trend between A(Mgb/(Fe)) and age 
identified in Fig. 1181 below. A remarkable conclusion from 
the comparison between the right- and left-hand panels in 
Fig. 1171 is that the chemical composition of early-type galax- 
ies appears to depend primarily on dynamical mass rather 
than stellar mass. 

It is also of interest to examine the relations between 
age, metallicity and a/Fe ratio. These are shown in Fig. 1181 
for our primary sample of early-type galaxies. In each panel, 
lines of different colours indicate the median trends fol- 
lowed by galaxies in different bins of velocity dispersion, 
from cry ~ 80 km s _1 (purple) to ~ 300 km s _1 (yel- 
low). Fig. 118b shows that, at fixed velocity dispersion, age 
and metallicity are stro ngly anticorrela t ed for galaxies older 
than log(t r /yr) = 9.8. iBernardi et alJ (120051) . who assume 
that the scatter about the CMR is caused entirely by age 
variations, suggest that young galaxies may be more metal- 
rich than older galaxies at fixed velocity dispersion. Other 
previous studies also invoked an age-metallicity anticorrela- 
tion to explain the tightne ss of the colour-magnitude and 
the Mg2-ov relations (e.g.. Ij0r gensenlll999t iFe rreras et alJ 
| l999l:lTrager et aljEoOOal : iTerlevich fc Forbesbool but see 
R^uMschlie^^t^ni200^1) ~^However. the slope of this anticor- 
relation in Fig. 118b is consistent with the slope of the age- 
metallicity degeneracy for individual galaxies quantified in 



Paper I using the same sample of early-type galaxies. Thus, 
we cannot conclude here whether this anticorrelation is real 
or induced by correlated errors, but we note that the ap- 
parent scatter in age of massive, old ellipticals is consistent 
with the measurement error and the scatter in metallicity is 
quite small. 

Fig. 118b shows the relation between A(Mgb/(Fe)) and 
metallicity for the galaxies in our sample. Interestingly, even 
though both the colour-magnitude and the Mg2-ov rela- 
tions imply an increase in both total metallicity and a/Fe 
with mass, we do not identify any significant correlation be- 
tween these two quantities (but see Fig. 119b below). This 
is true for the sample as a whole and for subsets of galax- 
ies with fixed velocity dispersion. Instead, A(Mgb/(Fe)) ap- 
pears to correlate more strongly with light-weighted age for 
log(t r /yr) > 9.7, independent of galaxy velocity dispersion 
(Fig. 118b ). Old galaxies appear to have larger a/Fe ratio 
than younger ones. This may be expected from the late en- 
richment in iron by Type la supernovae in galaxies which 
completed their star formation more recently, if the a/Fe 
ratio reflects the star formation timescale. We might be 
tempted to relate star formation efficiency (and hence a- 
enhancement) to stellar surface mass density, if the scaling 
between star formation rate and gas surface mass density in 
spiral galaxies ca n be extended to the progenitor s of early- 
type galaxies (e.g. ISchmidtJll959l : lKennicuttlll998l) . Fig. [lib 
shows A(Mgb/(Fe)) as a function of surface stellar mass 
density for the galaxies in our sample. The lines of different 
colours indicate the median trends followed by galaxies in 
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Figure 18. Correlations between stellar metallicity, light-weighted age and a/Fe-estimator A(Mgb/(Fe}). The solid lines represent the 
median relations for galaxies in different bins of velocity dispersion, centred on log((Ty /kms — 1 ) = 1.9,2,2.1,2.2,2.3,2.4,2.5 (increasing 
from purple to yellow). The grey scale indicates the density of galaxies with respect to the total at fixed abscissa. 



different age bins, from t r ~ 3 to 10 Gyr. We do not detect 
any significant correlation between A(Mgb/(Fe}) and sur- 
face stellar mass density. For completeness, in Fig. 119b . we 
show again A(Mgb/(Fe)) as a function of metallicity, as in 
Fig. 118b (the grey scale in this case indicates the distribution 
of galaxies as a function of A(Mgb/(Fe}) at fixed metallic- 
ity). There appears to be a mild correlation between a/Fe 
and metallicity for galaxies in fixed age bins, in agreement 
with the results of lTraeer et alJ j2000bl) . 



4 SUMMARY AND CONCLUSIONS 

We have analysed a sample of 26,003 early-type galaxies se- 
lected from the SDSS DR2 on the basis of their light concen- 
tration. Light-weighted ages, stellar metallicities and stellar 
masses for this sample were previously derived through the 
comparison of a set of carefully selected spectral absorp- 
tion features with a comprehensive library of high-resolution 
population synthesis models, encompassing the full range of 
physically plausible star formation histories (paper I). In ad- 
dition to these physical parameters, we have considered here 
an empirical estimate of the cv-elements-to-iron abundance 
ratio, given by the offset A(Mgb/(Fe)) between the observed 
ratio of the Mgb and (Fe) indices of a galaxy and that pre- 
dicted by the best fitting model in the library (which has 
scaled-solar abundance ratios). We have used these data to 
investigate the physical origin of two well-known scaling re- 
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Figure 19. a/Fe-estimator A(Mgb/(Fe)) plotted against (a) stel- 
lar metallicity and (b) surface stellar mass density within the r- 
band Petrosian half-light radius Rso.r- The solid lines represent 
the median relations for galaxies in different bins of light- weighted 
age, centred on log(t r /yr) = 9.5, 9.6, 9.7, 9.8, 9.9, 10, 10.1 (increas- 
ing from purple to yellow). 



lations for early-type galaxies: the colour-magnitude and the 
Mg2-ov relations. 

Our analysis demonstrates unambiguously and with 
unprecedentedly good statistics that both the colour- 
magnitude and the Mg2-ov relations are primarily se- 
quences in galaxy stellar mass. At increasing stellar mass, 
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as traced by either luminosity or velocity dispersion, the in- 
creasing colour and Mg- absorption line strength along the 
relations reflect an increase in both total metallicity and 
a/Fe ratio. Moreover, the galaxies in our sample cover a 
range in age of about 3 — 4 Gyrs, with more massive galaxies 
being on average older than low-mass galaxies. While at high 
masses early-type galaxies have the same mean age (and a 
small scatter in metallicity), at lower masses there is an in- 
creasing spread toward younger ages. This age spread domi- 
nates the scatter about the observed relations at low masses, 
in the sense that younger galaxies deviate toward bluer 
colours and lower Mg2 index strengths than older galaxies 
of the same mass. These results are consistent with the con- 
clusions from previ ous studies based on smal l er samples of 
early - type galaxies (iKoda ma fc Arimotdll 997UColless et alJ 
Il999l: IVazdekisll200ll : IWorthev fc Colloberdl2003h . In addi- 
tion, we find that the scatter in metallicity at fixed stellar 
mass contributes significantly to the scatter about the two 
observational scaling relations, in particular at high masses. 

We have checked that our main conclusions are not 
affected by possible dust effects and that they are robust 
against sample selection. In particular, the possible con- 
tamination of our sample by bulge-dominated star-forming 
galaxies, which could amount to ~10 percent, does not sub- 
stantially affect on our results. 

For a small subsample of 1765 galaxies we used 
information on enviro nmental density available from 
iKauffmann et all l)2004h to explore the dependence of the 
observed scaling relations and the stellar physical parame- 
ters on environment. We have found a small but detectable 
difference in the zero-point of the two relations, in the 
sense that early-type galaxies in dense environments tend 
to have redder colours and stronger Mg2 absorption indices 
than galaxies in low-density regions, at fixed luminosity or 
velocity dispersion. We also find a systematic increase in 
the scatter about both relations from high to low densi- 
ties. These variations appear to be induced by small differ- 
ences in the light-weighted age and metallicity of galaxies 
located in different environments. While galaxies with sim- 
ilar mass have the same element abundance ratio regard- 
less of environment, there is an increasing spread toward 
younger ages and lower metallicities in low-density environ- 
ments. At fixed stellar mass, early-type galaxies in dense 
regions are on average 0.02 dex older and more metal-rich 
than early-type galaxies in low-density regions. We note that 
these trends are very small and it will be worth re-examining 
them when a better statistics is available. If confirmed, these 
results are in agreement with previous studies indicating 
that early-type galaxies in clusters started to form stars ear- 
lier than, but on the same timescale as early- type galaxies 
in the field (e.g. IThomas et alj|200l iBernardi et"al]|2006t 
IClemens et alJl2006>) . 

We have also studied the dependence of the stellar mass 
on the dynamical mass estimated within the r-band Pet- 
rosian half-light radius of a galaxy. The relation is well de- 
scribed by a power law of exponent 0.783 ± 0.019 (equa- 
tion^IJl, implying a decrease in the stellar-to-dynamical mass 
ratio from low- to high-mass galaxies. The correlations of 
physical parameters with the dynamical mass estimated in 
this way suggests that metallicity and element abundance 
ratios in early-type galaxies are more fundamentally related 
to dynamical mass than to stellar mass. 



The increase in total metallicity with dynamical mass 
favours the classical interpretation of the colour-magnitude 
and Mg2-oy relations in terms of supernova -driven winds 
fe.g. iTraeer et ablEoOOal IThomas et al.ll2005l) . To account 
for the simultaneous increase in total metallicity and a/Fe 
ratio with mass, galactic winds should occur at early times, 
i.e. prior to the onset of Type la supernovae (on a timescale 
of a few Gyr), which are the main contributors to Fe-peak 
elements. The winds will be more effective in removing a ele- 
ments (produced by Type II supernovae on a time scale of ~ 
10 s yr) from low-mass galaxies with shallow potential wells, 
while they should not reduce significant ly the fraction of pri- 
mord ial gas and hence star formation iMac Low fc Ferraral 
1999). Much observational evidence has been accumulated 
for the importance of galactic outflows in galaxies with 
masses up to at least 10 10 Mp, jLehner t fc Heckmanlll996l : 
iHeckman et al.ll2000l IPettini et alJl2000h . Our results, if in- 
terpreted in terms of galactic winds, indicate that even more 
massive galaxies (with masses up to ~ lO n M0 in stars) 
have been affected by the ejection of metals through galac- 
tic winds. 

The above scenario, however, cannot account alone for 
the observed values of a/Fe in early- type galaxies. In a 
galactic- wind scenario, massive galaxies are predicted to 
have solar a/Fe abundance ratios, while low-mass galaxies, 
which lose a elements at early times, should have lower than 
solar a/Fe ratios. Instead, galaxies with stellar masses less 
than about 3 x 1O 1O M0 (corresponding to velocity disper- 
sions less than ~ 100 km s -1 ) are observed to have nearly 
solar a/Fe ratios, while this ratio increases to super-solar 
values in more massive galaxies (the quantity A(Mgb/(Fe)) 
reaching values around 0.3 in galaxies with stellar masses 
near 3 x lO n M ; see Fig.tlTl. 

The super-solar a/Fe abundance ratios of massive early- 
type galaxies suggest that these formed on a relatively short 
timescale and/or have/have had an IMF skewed towards 
high-mass stars. An IMF enriched in massive stars will pro- 
duce a larger ratio of Type II to Type la sup e rnova e, and 
hence a larger a/Fe ratio. iNagashima et ail (|2Q05) have 
shown that a top-heavy IMF during the burst ignited by 
the major merger that formed an elliptical galaxy can re- 
produce the observed range in a/Fe ratios. However, none 
of the models they explore yields the observed correlation 
of a/Fe with velocity dispersion (a model in which thermal 
conduction prevents the gas from cooling at the centres of 
massive halos is able to produce an increase in the a-element 
abundance, but not the a/Fe ratio, with mass). 

An interpretation of the a/Fe ratio in terms of the 
star formation timescale is supported by the correlation 
we find between A(Mgb/(Fe)) and light-weighted age, in- 
dependent of mass (Fig. I18H . This suggests that galaxies 
with longer star formation timescales (and thus with more 
recently formed stellar populations) have lower a/Fe ratio 
than galaxies formed on shorter timescales, because they 
had time to recycle the Fe-peak elements ejected by Type la 
supernovae. We also find that light- weighted age increases 
with stellar mass with negligible scatter at masses above 
10 1 Mq (Fig. I17t 0. Early- type galaxies less massive than 
about ~ lO n M0 display an extended tail toward younger 
ages, the mean age declining markedly with decreasing mass. 
This suggests either that low-mass galaxies formed more re- 
cently than high-mass galaxies, or that they have a more 
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extended star formation history (consistent with their solar 
a/Fe ratios). 

Our results represent further evidence for a shift in stel- 
lar growth toward less massive galaxies in recent epochs 
JCowie et all Il996t [De Lucia et all 12004 iKodama fc et alJ 
l2004l:lYi et all2005trTreu et all2005fl . This 'downsizing' sce- 
nario may appear at odds with the expectations of origi- 
nal hierarchical models of galaxy formation. However, ob- 
servations and the hierarchical paradigm can be recon- 
ciled if detailed physics of feedback from supernovae, ac- 
tive galactic nuclei o r thermal conduction i s introduced (e.g. 
[ Benson et all l200.lt lOranato et"aH 12004 iNaeashima et alJ 
2005). These sources of feedback could inhibit star formation 
on timescales short enough for the bulk of the star forma- 
tion to be completed before Type la supernovae can sub- 
stantially increase th e iron abundance in massive galaxies. 
ISoringel et, al.l l|2005h have shown that, in major mergers 
of spiral galaxies hosting supermassive black holes, AGN 
feedback provides a mechanism that can quench star forma- 
tion on short timescales. This mechanism is more efficient 
in the most massive early-type galaxies and leaves dwarf 
spheroids almost unaffected. We also note that the short 
star formation timescales (i.e. high formation redshifts) of 
massive ear ly- type galaxies do no t preclude longer assembly 
timescales foe Lucia et aljl2006l) : massive early- type galax- 
ies could appear old even if they assembled relatively re- 
cently. In this context, the new constraints derived here on 
the physical origin of the colour-magnitude and Mg2-ov re- 
lations for early-type galaxies represent a valuable reference 
for future models. 
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